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THIN  CYLINDRICAL  SHELL  UNDER 
INTERNAL  PRESSURE  AND 
CONCENTRATED  NORUAL  LOAD 
J.  E.  Goldberg 


ABSTRACT 


A  Method  is  developed  in  the  preeent  report  lor  calculating  the  noraal  eis- 
nleceBenta  in  a  very  thin,  pressurised  cylindrical  shell  which  is  subjected 
to  a  concentrated  noraal  lead*  The  shell  is  treated  as  a  aeabrane  and  the 
solution  is  obtained  in  the  fora  of  a  single  or  one-diaenaional  Fourier 
series,  the  longitudinal  coordinate  being  taken  as  the  independent  variable 
of  the  trigonoaetrie  teras,  with  the  cireuaf erential  com  iinate  being  ab¬ 
sorbed  into  the  Fourier  coefficients. 

The  solution  given  in  the  present  report  is  an  alternative  to  a  solution 
given  in  a  previous  report  in  the  fora  of  a  double  Fourier  series* 

The  report  concludes  with  a  miner i cal  exaaple. 


INTRODUCTION 


The  present  report  concerned  with  the  problea  of  determining  the  dis- 
placeaents  in  a  very  thin,  pressurised  cylindrical  shell  which  is  subjected 
to  an  "alaost  concentrated"  noraal  load. 

In  a  previous  report*,  a  aethod  was  given  for  deteraining  the  displaceaents 
under  the  loading  stated  above  by  first  expanding  the  "concentrated"  lead 
into  a  double  Fourier  series.  It  was  assuaed  that  the  displaceaents  could 
also  be  represented  by  a  double  Fourier  series  of  the  saae  trigonoaetrie 
functions,  and  it  was  shown  that  a  definite  relation  existed  between  the  co¬ 
efficients  of  the  loading  series  and  the  coefficients  of  J i splaceaent 

series. 

Unfortunately,  the  series  n.v~;ved  in  the  aethod  previously  described  is 
likely  to  converge  quite  slowly  in  certain  actual  cases  and,  in  these  i-ses, 
it  becoaea  necessary  to  coapute  a  large  nuaber  of  teras  in  order  to  obtain 
the  required  accuracy. 


METHOD  OF  /•J.'A-.fSIS 


We  present  now  a  aethod  of  analysis  which  could  reduce  rater xa\ ly  the  aaount 
of  coaputation  involved  in  obtaining  a  solution  to  the  ;>rableu  stated  above. 


1.  Goldberg,  J.  E. ,  "Theory  of  Meabranes,  Deflections  *  d  Jtr«:ss*s  in  Mem¬ 
branes  unde”  Concentrated  Noraal  Load",  KRR-AN-05F,  ’5  July  *90! 
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We  Miuat  that  the  oorail  load  P  ia  located  on  the  surface  of  the  shell  at 
(in  cylindrical  coordinates)  ©  ■  0,  *  -  e  where  e  is  aeasured  froa  one  end. 
Instead  of  a  single  load  P  at  that  location,  wo  shall  assume  that  we  have, 
at  x  *  c,  two  loading  cases  as  follows: 

(a)  Symmetrical  loads  each  of  magnitude  ?/ 2  located  at  ©  =  0 .IT 

(b)  Anti-symmetrical  loads  earh  of  aagnitude  P/2  located  at  &  -  0,T T 

Clearly,  the  loading  obtained  by  superimposing  loadings  (a)  and  (b)  is 
equivalent  to  the  original  single  load  P  at  d  .  0. 


We  will  now  imagine  that  the  shell  is  divided  into  two  aeai-cyl indri cal 
shells  by  cutting  at  6  ■  0  and  0  *1r  ,  and  we  assume  that  the  several  loads 
P/2  which  have  been  placed  on  the  lines  at  0  «  0  and  ©  m'TT  have  also  been 
split  into  two  equal  components. 


drieal  "membrane"  which,  in  addition  to  an  internal  pressure  is 
along  its  straight  edges  as  shown. 


loaded 


We  as  "ise  that  the  internal  pressure  ia  such  as  tc  cause  s  hoop  membrane 

forev  (i.e.,  N  -  pR)  and  an  axial  membrane  force  S  /2.  We  shall  take 

*  o  o 

w  to  he  the  additional  deflection  (i.e.,  in  addition  to  that  due  to  inter¬ 
nal  pressure  alone)  due  to  the  applied  load  P. 

Within  the  boundaries  of  each  semi-cylinder  tbe  deflection  oust  satisfy 
the  equation 


ik  ii- ♦  -  o 

R2  3G2  2  3x2 


(1) 


We  will  assume  that  the  deflection  may  be  taken  in  the  form  of  a  series 

w (x,6)  -  2  W  ain  nlfcc/L  (2) 

n*  i  n 

where  the  factors  W  -  W  (6),  that  is,  the  coefficients,  are  functions 
•f  the  angular  position?  This  series  will  give  zero  normal  displacement  at 
the  circular  edges  of  the  aeai-cyl indrical  shell. 

We  will  consider  the  implications  of  one  term  of  tbe  series  in  Equation  (2), 
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that  is,  we  will  assuae  that 


w  a  w  -  If  sin  n'Tr  x/l. 
n  n 

Substitution  of  Equation  (3)  into  Equation  (1)  would  lead  ts 

o  2 

d“W  -  k  W  *  0 
n  n  n 

O 

d~“ 

C  2  2 

where  k  -  R  ,n TT  \ 
n 


The  general  solution  of  Equation  {■♦)  is 


V  -A  sinh  k  6  +  H  cosh  k  0 
n  n  n  n  it 

Along  the  edge  &  »  0, 

»  A  k  sin  nTI* 
n  n  n  _ 

R  P©  R  L 


(3) 

(4) 

(5) 

(«) 


And  therefore  the  distributed  noraal  (i.e.,  radial)  force  along  this 
edge  is 


f  »  S  A  k  sin  riTfx 
n  o  n  n  —r— 

R  L 


(7) 


Let  us  now  consider  a  "concentrated"  load  •/  -  YjA  which  is  applied  at  the 
point  x  =  c  on  the  edge  6>  0.  This  "co*  centrated"  load  nay  be  assuaed  to 
be  distributed  over  a  saall  length  LV  in  the  direction  of  the  edge,  and  we 
nay  therefore  say  that  Q  is  e«(u i  v- 1  e>. 4.  !c  a  distribute*!  loading  along  the 
edge  defined  as  follows: 

SO  0<x£(c  -  e) 

Q/2 «  (c  -  e)£x^(c  «•  r)  (8) 

0  (c  ♦  e)fx£L 

Such  a  loading  nay  be  represented  by  a  Fourier  sertea: 

OO 

Q —  9  *,n  nTTx/L  (9) 

n-1 

'A'e  have  now  to  determine  the  coefficients  qn  by  a  standard  technique.  t»e 
substitute  Equation  (8)  into  Equation  (9)t  multiply  both  sides  of  the  re¬ 
sulting  equation  by  sin  iTTx/L  dx  (where  i  is  an  arbitrary  integer)  and  in¬ 
tegrate-  -e  observe  that  substitution  of  Equation  (8)  into  the  left  side 
of  Equation  (9)  and  performing  the  operations  mentioned  above  would  yield 

'J  q  sin  iTTx  dx  *  0  ♦ 

°  L  2e 

c-e 

If  e  is  taken  sufficiently  small,  the  integrand  of  the  term  on  the  right 
side  remains  essentially  constant^,  and  we  obtain 


C  +  e 


sin  iTTx  dx 
L 
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X 


L  pc+e 

q  sin  i1Tx  dx  -  I  -in  dx 

.  L  2*4-«  1 

L 

q  aia  lflx  dx  »  Q  ain  i*iTc 
L  L 


(10) 


Perforaing  the  saae  operationa  upon  the  right  aide  of  Equation  (9)  we  have, 
by  uae  of  the  orthogonality  relatione 
\L 

ain 

*  L  L  iJ‘  L 


fcj  iri.wuu. 

P (2^  ***»  n]Tx)  sin  iTtx  dx  ■  Pq.  sin**  lflx  dx 

J,  L  L  Jo  1  L 


-  OjL 
2 

Equating  (10)  and  (ll)  and  replacing  i  by  n  yielda 

*  2£  sin  n*Rc 
L  L 


(n) 


(12) 


Hence,  the  concentrated  load  y  («P/4)  aay  be  repreaented  by  the  aeries 

oo 

(13) 


Q-~'2«l/L  ^  (sin  riHc^in  riTTxj 


Thus  a  concentrated  load  correaponds  to  the  auperpoaition  of  a  aeries  of 
sinusiodally  varying  distributed  loads. 

Ve  aay  now  eoapare  one  tera  of  the  series  in  Equation  (13)  with  the  sinu¬ 
soidal  tera  of  Equation  (7)  and  we  see  that  the  two  teras  are  identical  if 


A  »  -  2 OR  ain  riffc 

“  LN  k  L 

o  n 

Recalling  that  Q  -  P/d,  and  using  Equation  (i)  this  becones 

A  « _  -F  sin  rifle  (id) 

n  ^2HN  n  L 

o 


Jr  r 

sin  itfx  dx  »  QL  Icos  i*TT(c-e }  -  cos  i1T(c  *  fe) 
c-e  L  2ei7T  L  L  L 

*  QL  sin  iTfc  sin  iffe  ^  QL  /sin  iTTcv  iHo  if  iflTe/L  is 

eiTT  L  L  eilT  '  L  '  L 

fairly  saall,  say  one-tenth  of  a  radian  or  less.  lienee,  if  we  can  "cut 
off"  our  series  at  a  point  where  this  restriction  is  satisfied,  the  re¬ 
sults  should  be  acceptable. 
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We  return  no*  to  consideration  of  Equation  (6),  the  general  solution  of  the 
ordinary  differential  equation  for  W  .  If  the  loads  P/4  are  applied  ijfue- 
trically,  then  the  deflections  alongC the  edges  &  -  0  and  0  *  TT  are  equal. 
Hence,  for  this  ease 

B  «  A  sinh  k  TT  ♦  B  cosh  k  T r 
n  n  n  n  n 

or  B  •  A  sinh  k  T"’  /( 1  -  cosh  k  Tl)  (symmetrical)  (15) 


If,  on  the  other  hand,  the  loads  p/4  are  applied  antisyaaetrieal ly,  the  de¬ 
flections  along  these  edges  are  equal  in  magnitude  but  opposite  in  sign. 
Hence  for  the  antisyaa »trieal  case 


or 


B  m  -A  sinh  k  /TT  -B  cosh  k  TT 
n  n  n  n  n 

B  ■  -A  sinh  k  Tr  /( 1  +  cosh  k  Tf)  (anti  syaaetr  ical ) 
n  n  c  n 


(16) 


A  concentrated  load  P/2  at  ©  »  0,  x  »  c  on  the  seai-shell  (corresponding  to 
a  single  load  P  on  the  coaplete  shell)  was  shown  to  be  equivalent  to  symae- 
trical  and  antisyaaetrical  loads  P/4  on  the  seai-shell.  Therefore  we  eoa- 
bine  both  Equationr  (15)  and  (16),  substitute  into  Equation  (6).  thence  into 
Equation  (3),  and  upon  using  Equation  ( 1 4 )  and  Equation  (2)  we  obtain  the 
following  expression  for  the  additional  noraal  deflections  of  a  very  thin 
pressurized  cylindrical  shell  when  subjected  to  a  concentrated  noraal  load 
P  applied  at  Q  *  0,  x  *  c: 


*(*,e) 


-  47P 
2T1N 


oo 


n 


sioaiix) 

’I T 


or 


w(x,G) 


-  y 

o  n-1 


1  cosh  k  (TT -  0) 

— _ n _ _ 

n  sinh  k  tt 
n 


sin  nfTc  sin  nTTx 
L  ~L 


(17) 


where 


b  ir 

n 


DISCUSSION 

The  coefficients  in  Equations  (17)  are  easily  evaluated  either  by  use  of  a 
table  of  hyperbolic  functions  or  by  use  of  the  fundaaentai  exponential  foras 
of  these  functions  and  a  table  of  powers  of  e.  When  the  argument  of  these 
hyperbolic  functions  becoaes  large,  asymptotic  foras  aay  be  useful.  For 
example:  sinh  0/eoeh  Q  -*■  1  -  2e“^e  as©  becoaes  wery  large. 

Equations  (17)  are  valid  when 

(a)  the  shell  is  very  thin 

(b)  the  internal  pressure  is  high 

(c)  the  ’’concentrated"  load  is  saall 
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When  these  restrictions  are  violated,  the  deflections  given  by  these  ecua- 
tions  generally  are  larger  than  the  actual  deflections.  Rnles  for  estab¬ 
lishing  the  range  of  validity  of  the  membrane  theory,  end  for  estimating 
the  secondary  or  additional  stresses  in  teras  of  the  displaceaente  given 
by  the  above  eqnations  were  presented  in  the  previous  report. 

We  will  not  apply  Equations  ( 1 7 )  to  the  point  directly  under  the  center  of 
the  lead  since  a  theoretically  concentrated  load  would  result  in  an  infinite 
deflection  at  the  lead  point.  Actually  however,  we  do  not  deal  with  theo¬ 
retically  concentrated  leads  but  with  loads  which  are  distributed  over  soae 
finite  area.  If  the  area  is  snail.  Equation*  (17;  should  give  reasonably 
goed  results  beyond  the  edge  of  the  actual  contact  area. 

It  would  be  well  to  avoid  sharp  edges  on  the  attachment  which  transits  the 
load  to  the  shell.  If  possible  the  slope  at  the  edge  of  the  fitting  should 
not  he  less  than  the  calculated  slope  of  the  shell  at  the  edge  ef  the  fit¬ 
ting.  Also,  the  curvature  of  the  fitting  should  not  exceed  that  which,  when 
applied  to  the  shell,  causes  extreae  fiber  stresses  of  a  magnitude  such  as 
to  raise  the  stress  level  above  the  design  allowable  stress. 


ILLUSTRATIVE  EXAMPLE 


As  an  illustrative  exaaple,  we  take  a  very  thin  cylindrical  shell  subjected 
to  a  positive  internal  pressure  ef  sufficient  magnitude  as  to  place  the 
shell  in  the  approximate  condition  of  a  aeabrane.  The  shell  is  subjected 
also  to  a  "concentrated"  normal  load  P  located  at  c  ■  L/2.  The  necessary 
data  include: 

R  -  40  inches  k  -  R'ff  n  -  .740480  n 

“  VsT  l 

L  «  120  inches  k  "rf  •  2.326287  □ 

n 

We  will  calculate  the  normal  displacement  at  the  point  x  -  L/2,  ©  -  0.1. 
Thus,  if  the  fitting  which  transmits  the  load  to  the  shell  has  a  dimension 
of  about  eight  inches  in  the  circumferential  direction,  is  suitably  shaped 
both  in  profile  and  plan  (probably  with  an  elliptical  planform  having  its 
major  a*jis  in  the  circumferential  direction),  and  is  itself  properly  de¬ 
signed,  the  caleulated  deflection  would  be  expected  te  pertain  to  the  edge 
ef  the  fitting. 
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te  observe  tnat  since  sin  c  "'7  c/L  »  0  for  all  eves  integral  values 
of  n,  it  is  necessary  to  calculate  terms  only  for  the  odd  integer*. 
Furthermore,  sines  c  *  L/2  and  x  *  L/2 ,  then  sin  nlT c  sin  nTT-  x  * 

L  L 

for  the  odd  integers.  S<  note  also  that  coth  k  TT  becomes  unity  (to  the 
sixth  decimal  place)  when  n  «  5  and  remains  at  ?h*>- value  for  larger  n. 
Consequently,  for  n  equal  to  or  larger  than  5, 

— k  B 

n  F  *  coth  k  iT  cosh  k  S  -  slab  k  £  *  e  n~ 

j  n  n  n 

Obviously  this  simplification  is  possible  only  when  k  is  sufficiently 
large.  Aside  from  this  simplification,  the  series  may  be  expected  to 
converge  nicely  if  6  (which  is  a  coordinate  of  the  point  where  the 
deflection  is  being  computed,  not  the  point  where  the  load  is  applied) 
is  not  too  small,  i.e.,  is  not  too  close  to  the  generator  upon  which 
ths  load  is  applied. 


The  calculations  are  carried  out  in  the  table  below.  Because  of  the 
simplification  mentioned  above,  sinh  k  9  and  cosh  k  9  are  cot 

D  fi 

tabulated  for  n  greater  than  5,  and  the  value  of  nil  for  these  values 
of  n  is  taken  directly  from  a  table  of  exponentials  as  9. 

The  last  column  of  the  table  should  be  multiplied  term  by  term  by  sis  a  TT  c 
sin  nTT  x  before  summing.  Since  this  product  was  shown  to  be  unity  L 
L 

(for  odd  values  of  n)  for  the  particular  loan  point  and  deflection  point 
under  consideration,  this  step  has  been  omitted.  In  the  general  case, 
obviously,  a  colon  should  be  provided  for  this  calculation. 


n 

k  TT 

c 

cth  k  TT 
n 

k  9 
n 

snh  k  © 
n 

esh  k  9 
n 

sl 

s 

» 

3 

i 

2.326287 

1.019258 

.074116 

.074116 

1.002742 

.947690 

.947390 

3 

6.978862 

1.000001 

.222144 

.223976 

1.024775 

.800800 

•266933 

5 

11,631436 

1.000000 

.370240 

. 37S757 

1.069325 

.69055® 

.133114 

7 

20.936585 

.518336 

.619035 

.033434 

9 

.666432 

.513533 

.057060 

11 

.814528 

.442348 

.040247 

13 

.962624 

.  381690 

.029376 

15 

1.110720 

.329322 

.021955 

17 

1.25SS16 

.233990 

.016705 

19 

1.405912 

.244399 

.012839 

21 

1.555008 

.211139 

.010057 

23 

1.703104 

.132113 

.007918 

25 

1.851200 

.15704® 

.006232 

27 

1.999296 

.135431 

.005016 

29 

2.147392 

.115733 

.004027 

1.552903 
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The  above  result  shows  that,  for  a  sufficiently  thin  shell  having 

the  diaenaions  stated  above  and  subjected  to  a  sufficiently  high  internal 

pressure,  the  noraal  deflection  at  the  point  x  »  L/2,  9  «  0*1  would  be 

*  ■  WTT-(1-6529) 


If  the  shell  and  loads  Jo  not  satisfy  the  stated  conditions,  a  solution 
of  the  above  type  aay  be  taken  as  a  first  approxiaation* 
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